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Two-leg spin ladders have a rich phase diagram if rung, diagonal and plaquette couplings are 
allowed for. Among the possible phases there are two Haldane-type spin liquid phases without local 
order parameter, which differ, however, in the topology of the short range valence bonds. We show 
that these phases can be distinguished numerically by two different string order parameters. We 
also point out that long range string- and dimer orders can coexist. 
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I. INTRODUCTION 

The properties of spin ladder systems attracted consid- 
erable interest in the past decade J' It was shown^ that 
even-leg ladder models develop a gap in their excitation 
spectrum if the rung coupling is ferromagnetic. This re- 
sult could be considered as a generalization of Haldane's 
conjucture^ about the existence of a finite gap in integer 
spin Heisenberg chains, since for strong enough ferromag- 
netic rung coupling the spin ladder with spin-S* on the 
legs is equivalent to a spin chain of spin-2S'. In this case 
the ground state is an AKLT like state^ in which short 
range valence bonds couple spins on neighboring rungs. 
Later it was demonstrated numerically,^ that a gap is 
generated for antiferromagnetic rung coupling as well, in 
which case the singlet valence bonds are formed predom- 
inantly on the rungs. It may be called a rung singlet 
state. Both are Haldane-type spin liquids in which there 
are no broken local symmetries, and the excitations are 
coherent magnons with a finite gap. 

Similar situation occurs if the two legs of the ladder 
are coupled not by rung couplings, but diagonally or 
by plaquette couplingsn^ When these couplings compete 
with each other, non-Haldane spin liquid phases, char- 
acterized by the absence of coherent magnon excitations 
and spontaneous dimerization,M.may occur due to frus- 
tration effects. Although the thermodynamic properties 
of Haldane and non-Haldane spin liquids are identical, 
the correlation functions differ due to the difference in 
the excitations spectrum. The two phases could also be 
distinguished by a local dimer order parameter. 

Beside these spin liquid phases several other phases 
may occur, like ferromagnetic or incommensurate phases. 
When drawing the phase diagram, it was thought for 
some time that the concept of Haldane-like spin liquids is 
unambiguous for spin ladders, and therefore the AKLT- 
like and rung singlet phases are identical not only thermo- 
dynamically, but also in the sense, that in a large param- 
eter space, where rung, diagonal and plaquette couplings 
are all taken into account, one could go from one to the 
other without any discontinuity. In fact, it was arguec^ 
that there is a controlled way of changing the parame- 
ters gradually and interpolate continuously between two 
limiting cases belonging trivially to one or to the other 
phase. The procedure contains a shift of one of the legs 



of the ladder by one unit, thereby transforming a rung 
singlet into a diagonally situated valence bond, charac- 
teristic of the the AKLT-like state, or vice versa. How- 
ever, as was pointed out by Kim et al^ this shift changes 
the topology of the system drastically. When we start 
from an open ladder, the above transformation deforms 
the usual straight-end ladder into a "ladder" whose ends 
are cut diagonally. Such a deformation has serious con- 
sequences on the degeneracy of the ground state for open 
boundary condition. In a usual ladder the AKLT state 
has a fourfold degeneracy due to the unpaired end spins, 
while the rung singlet state is unique. On the other hand, 
in the diagonally cut ladder the ground state is unique 
for the AKLT state, and fourfold degenerate in the rung 
singlet state. This is best seen in Fig.^a) and (b), where 
the typical valence-bond configurations of the AKLT-like 
state and the rung singlet state are shown. In Fig. ^c) 
we also show a typical valence bond configuration of the 
dimer phase. In this phase the invariance of translation 
by one rung is broken spontaneously. 

While in the AKLT state the number of valence bonds 
crossing an arbitrary vertical line between rungs is al- 
ways odd, in the rung singlet state it is always even. 
This is true not only in this special configuration but 
remains valid generally, as long as the range of valence 
bonds is short compared to the size of the system.^" The 
two states can be distinguished by a topological quan- 
tum number. Thus for topological reasons in a spin lad- 
der the two Haldane-type spin liquid phases are different 
and in the space of couplings they must be separated by a 
phase transition hyperplanef2*ii although extended gap- 
less (critical) regions between them cannot be excluded 
a priori either. 

Since the topological quantum number is not accessible 
to direct computation, a generalization of the hidden or- 
der parameter of the spin-1 chainiS was proposed,^ that 
could distinguish between the even and odd topological 
sectors of the Haldane phase. Although the relationship 
between the two new string order parameters and the two 
types of Haldane phases was made plausible by looking 
at special configurations, it was not possible to calculate 
their value, since the bosonization procedure applied to 
obtain the phase transition lines did not distinguish be- 
tween the two parameters. 

The aim of this paper is to do a numerical study of 
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the relationship between the possible spin liquid phases 
of the spin ladders and the string order parameter, us- 
ing the density matrix renormalization group (DMRG) 
procedureJ^ The setup of the paper is as follows. In Sec. 
II we give a short description of the ladder models and a 
brief account of what is known about the possible phases. 
The results of our numerical calculations are presented in 
Sec. III. Finally Sec. IV contains a brief summary. 



II. GAPPED PHASES IN LADDER MODELS 

Following the notations of Ref. IIVI we will write the 
most general form of the Hamiltonian for isotropic two- 
leg spin ladder models if only spins on the same and 
neighboring rungs interact, as 



spin ladders by defining the order parameter of AKLT- 
like Haldane state as 
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where 



hi^i+l — J\ Si ■ Ui+l + J{ Ti ■ Ti+l 

1 ^ ^ 1 , ^ 

+ Si ■ fi+i + J^fi ■ Si+i 

+Vii {Si ■ Si+i){n ■ n+i) 

+ Vrr [Si ■ Ti){Si+i ■ fi+i) 
+ Vdd (Si ■ Ti+i){Ti ■ Si+i) . 
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In Eq. (0) cTi and n are spin- 1/2 operators correspond- 
ing to spins sitting on the two legs of the ladder. The 
schematic plot of the spin couplings in ft-^.i+i between 
the spins on the two legs of the ladder are shown in Fig. 

m 

The terms with Ji and J( couple neighbouring spins 
on the same legs, the terms with Ji- and J^. spins on the 
same rung, while the terms with Jd and spins that are 
situated diagonally. In addition to these bilinear terms 
there could be three biquadratic terms with couplings Vn, 
and Vdd- 

In order to decrease the number of parameters from 
now on we only consider the case when J^i = Jr — J^, 
Ji, = J\ = J[ and Jd = Jd = ^he energy scale will 
be set by fixing Jl = 1. 

It is known from previous studies^ that all the interleg 
couplings Jr, Jd, . • . , Vdd, when acting alone, generate a 
Haldane-like phase. This is an AKLT-like phase with 
odd topology of the number of valence bonds if Jr < 0, 
Jd > 0, V^i > 0, < or Vdd > 0, while for opposite 
sign of the couplings a rung singlet like state with even 
topology is generated. 

As mentioned before, the goal of the present paper 
is to identify an appropriate order parameter that can 
characterize the spin liquid phases. Starting from the 
string order parameter of the S' = 1 chain proposed by 
den Nijs and Rommelse— , it is natural to generalize it to 
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where S*" is the total spin on rung i, Sf = af + t", 
a = x,y, z. 

Since in the extreme limits the pure rung singlet state 
can be related to the AKLT state by a shift of one of 
the legs, the string order parameter of the Haldane state 
with even topological quantum number is expected to be 
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where Sf is now the sum of two diagonally situated spins, 

ca a I a 9.15 

•^i — ^ ^i+l'^ 

In what follows we will calculate the order parameters 
for the z component only, they will be denoted as Oodd, 
and Oevon, respectively. We will show, that these two or- 
der parameters are mutually exclusive. Oodd is finite for 
the AKLT-like Haldane state, while Oovcn is finite for the 
rung singlet like Haldane phase. Moreover the order pa- 
rameter goes to zero continuously at the transition point 
when the transition is of second order, while it vanishes 
with a jump in case of first order transitions. 



III. NUMERICAL RESULTS 

To investigate this problem we performed numerical 
calculations applying the DMRG method^^ on the model 
defined by the Hamiltonian in Eq. Since the DMRG 
procedure is more accurate for systems with free ends, 
we consider our ladder model with open boundary condi- 
tion. In order to eliminate the fourfold degeneracy of the 
ground state in the conventional Haldane phase due to 
the free S* = 1/2 spins at the ends of the ladder, we ex- 
ploited the left-right reflection symmetry of the system. 
Out of the four degenerate ground state we have selected 
the one with S^^^ = and positive parity. The calcula- 
tions were done for chains with up to 150 sites keeping 
100 to 200 states per block. The truncation error varied 
in the range 10~^ to 10^^. In most of the calculations we 
used the finite-lattice method with two or three iteration 
cycles. 

In order to reduce the undesirable boundary effects 
we fixed site i in Eqs. Q and Q at the center of the 
chain and let site j run through half of the chain. The 
finite size scaling procedure was used to determine the 
thermodynamic limit of the quantitities measured. 

In what follows we present our numerical results for the 
string order parameters for various choices of the coupling 
constants. Since our aim was not to explore the full phase 
diagram, but rather to check how the string order param- 
eters Oeven and Oodd Vary as we go from one Haldane-like 
phase to another, topologically different gapped phase. 
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we calculated the string order parameters along specially 
chosen paths connecting points in the parameter space, 
where the model is known to belong to different phases. 

It is worth pointing out that since the ground state 
wave function and therefore the ground state expectation 
values of operators can be calculated more accurately 
than those of excited states, the order parameters provide 
a better way to obtain phase diagrams than investigating 
the energy spectrum. 



A. Haldane phases due to rung, diagonal or 
plaquette couplings alone 

The effect of Jr alone was analysed numerically for 
short systems in Ref. IIVI It was found there that for 
Jr < the order parameter Oodd is finite, while for Jr > 
the other order parameter Ocvon is finite. Our results 
with the DMRG method confirmed these results. 

From earlier works- it is known, that not only the rung 
but also the diagonal and plaquette inter-leg couplings 
are relevant perturbations which generate a gap and drive 
the system into phases with topological order. Depend- 
ing on the sign of the perturbation the emerging phases 
belong to one of the two possible topological sectors. 

We have therefore investigated the effect of these cou- 
plings by doing calculations at two specialy chosen values. 
The finite size scaling analysis of the results obtained for 
couplings Jd = ±1 or Jd = il indicated that irrespective 
whether only one or both diagonal couplings are present, 
Codd / and Ocvcn = for positive Jd or Jd, while we 
have Ocvcn 7^ and Oodd = for negative Jd or Jd. 

Similarly, calculations for the values ±1 of the cou- 
plings Va, Vrr and Vdd (the other couplings set to zero) 
convinced us that the even order parameter is finite, 
Ocvon > 0, for Vii < 0, or V„ > or Vdd < 0, while 
the odd order parameter is finite, Oodd > 0, for opposite 
signs of the couphngs, i.e., for Vii > 0, or V„ < or 
Vm > 0. 

Our results obtained for the string order parameters 
allows us to conclude that the topology of the valence 
bond structure can be identified directly by calculating 
the order parameters Oodd or Oovcn , since they mutually 
exclude each other. 



B. The Jr - Jd phase diagram 

Next we consider the subspace spanned by Jr and Jd . 
According to the recent bosonization study^-^ the phase 
transition line separating the two topological sectors is 
of first order in the quadrant Jr > 0, Jd > 0, while of 
second order in the quadrant Jr < 0, Jd < 0. In order 
to check the behavior of the string order parameter, we 
have calculated them along two paths in the (Jr, Jd) 
parameter space. 

Path 1 is defined by Jr = 1 — Jd with < Jd < 1, 
and path 2 is parametrized by Jr = — 1 — Jd with — 1 < 



Jd < 0. They are indicated by sohd lines on Fig.|3| 

Our numerical results for Oeven and Oodd obtained in 
the thermodynamic limit along path 1 is plotted on Fig. 
0] The lines connecting the points are only guides to the 
eye. It is apparent from the figure that Oovcn is finite in 
one particular region while Oodd in the other. 

As one sees from the figure the string order param- 
eters, extrapolated to the infinite ladder limit, have a 
jump discontinuity at the transition point. In order to 
further support the first order nature of the phase transi- 
tion we have also investigated the low lying energy spec- 
trum. We found that there is a sudden change in the 
degeneracy of the ground state at the transition point. 
In the case of open ladders the fourfold degeneracy of 
the AKLT phase associated with spin-1/2 end-spins dis- 
appear, and the ground state becomes unique as a rung- 
dimer-like structure emerges. Moreover, by examining 
higher lying energy levels we found that the energy gap 
remains finite throughout the transition. The ground 
state energy shows a cusp-like singularity, a characteris- 
tic feature of first order transitions, as is shown on the 
inset of Fig.^ Note, however, that for finite chains there 
is no observable level crossing in the ground state due to 
hybridization. 

Next we consider the behavior of the order parame- 
ters calculated along path 2 in the negative Jr and Jd 
quadrant. Our results are plotted in Fig. |S1 In this case 
again, only one of the order parameters is finite in a par- 
ticular phase, however, the gap vanishes continuously at 
the phase boundary, which is a characteristic features of 
second order phase transitions. It is worth noting that 
there is an order of magnitude difference in the scale of 
Fig. ^ and Fig. O as a function of Jd ■ 

Since the phase transition points were located on two 
special trajectories only, the phase boundary that sepa- 
rates the two types of Haldane-like phases is drawn by 
dashed line in Fig. to indicate that its exact shape was 
not determined. 



C. The role of plaquette couplings 

In order to consider the infiuence of the plaquette cou- 
plings we have chosen the special values K = Jb. = Vi\ = 
Vdd) Vrr = 0, since the phase diagram has already been 
studied for this parametrization.^ This choice has been 
motivated by the possible mapping of the ladder model to 
the spin-1 bilinear-biquadratic chain for K = 1.^ Here we 
calculate string order parameters along four lines going 
through the expected phase boundaries shown by solid 
lines in Fig. 

Calculations along path 1 have revealed that Oodd is 
finite in the AKLT-like Haldane phase, but it vanishes in 
the dimerized phase. On the other hand, Oeven becomes 
finite in the dimerized phase, thus the two phases are 
well distinguished by the order parameters. Examining 
the trajectory (path 2) going through the phase bound- 
ary between the rung singlet like Haldane and dimerized 
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phases, we have found that Ocvcn remains finite and Oodd 
scales to zero everywhere. This shows that the string 
order parameter does not uniquely characterize the Hal- 
dane phase, it is rather an indicator of the topology of 
the valence bond structure, since both the rung singlet 
and the dimerized phases have even topological quantum 
numbers. 

Thus when Ocvcn 7^ 0, the string order is not sufficient 
to classify unambiguously the ground state. Another or- 
der parameter, the so-called dimer order parameter has 
to be measured. Taking the difference of the local energy 
on neighboring bonds in the middle of the ladder this 
order parameter is defined as 

In the Haldane-type spin liquids D vanishes, but there 
exists an extended region in the parameter space, where 
the dimer order parameter scales to finite value in the 
thermodynamic limit, signaling a dimerized phase with 
spontaneously broken translational symmetry. 

Next let us consider the K > Q half-plane focusing 
on the Jd = half-line and its vicinity. Calculation 
of the order parameters along path 3, defined hy K ^ 
1 -I- Jd with — 1 < Jd < 1 revealed that Oeven is finite for 
Jd < 0, while Oodd is finite for Jd > 0. At Jd = all 
order parameters including D vanish as is shown in Fig. 
Since the string order parameters become finite for 
arbitrary small finite value of Jd the two types of spin 
liquids are separated by a critical line and the extended 
critical phase reported earliest does not exist. 

Finally, it can be seen on the figure that for Jd > 0.6 
both order parameters vanish and a new phase appears. 
In order to analyze the question whether the new phase 
develops just above the VBS point, we investigated the 
vicinity of the VBS point along path 4. We found that 
Oodd remains finite up to if 1.9 while Oeven vanishes 
for the whole regime. Investigating the low lying energy 
spectra we found a multiply degenerate ground state. In 
order to decide whether this phase extends to the Jd axis 
we made test calculations at Jd — 0.05, K = 1.5,2 and 
found a finite Oodd but vanishing Oeven order parameters. 
This suggests the schematic phase boundary shown by a 
dashed line in the figure. 

IV. CONCLUSIONS 

In the present paper we have studied the relationship 
between various gapped Haldane and non-Haldane spin 



liquid phases and the two different string order param- 
eters in the general two-leg spin ladder model, contain- 
ing rung, diagonal and plaquette couplings as well. The 
string order parameters have been calculated numerically 
using the DMRG method. 

It has been found that there are two kinds of Haldane- 
type massive spin liquids characterized by finite Ocvcn or 
Oodd string order parameters, respectively. It was shown 
that Oodd is finite in the AKLT-like phase, while Ocvcn is 
finite in the rung singlet like phase. Examining the tran- 
sition from the AKLT-like to the rung-singlet-like phases 
the type of phase transition can be first or second or- 
der depending on the path in parameter space, and this 
shows up in the abrupt or continuous vanishing of the 
order parameter. 

When plaquette couplings are allowed for, non- 
Haldane-likc spin liquid phases can occur, and the string 
order parameters are not sufficient to characterize the 
phases. Ocvcn remains finite in the dimerized phases, 
too. Thus the string order parameter reflects the topol- 
ogy of the valence bond structure and does not charac- 
terize uniquely the Haldane-type phases. A dimer order 
parameter can, however, be defined, which is finite in a 
non-Haldane-type dimer phases but vanishes in Haldane- 
type spin liquids. 

It was found previously that there exists a small range 
of the parameters where the combination of perturbing 
operators is irrelevant and the spectrum remains mass- 
less. Investigating the vicinity of Jd = 0, if > axis we 
have found that both string order parameters vanish on 
this critical line, but Ocvcn or Oodd becomes finite im- 
mediately as we move away from this axis. Therefore, 
we conclude that the two kind of spin liquids are sepa- 
rated by a critical line and not by an extended gapless 
regime reported earlier. In addition, around the Jd = 1 
line for K > l.i we have found a new phase. This phase 
is characterized by multiply degenerate ground state and 
vanishing order parameters. The detailed analysis of this 
region is left for future investigation. 
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Fig. 1 Path et al. 



FIG. 1: Representative valence-bond configurations without quantum fluctuations in diflFerent gapped phases of spin ladders, 
(a) The AKLT-like state, (b) rung singlet state, (c) dimerized state. State (a) is odd, while states (b) and (c) are even under 
our topological classification. 
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FIG. 3: The Jr-Jd phase diagram. Phase boundary is indicated by dashed line and the two trajectories of the calculation by 
sohd hnes. 
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Fig. 4 Path et al. 



FIG. 4: The N oo extrapolated values of the string order parameters Oovcn and Oodd calculated along path 1. The inset of 
the figure shows the cusp singularity of the ground state energy indicating a first order transition. 




Fig. 5 Fatli et al. 



FIG. 5: The Oeven and Oodd string order parameters calculated along path 2. 
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Fig, 6 Path et, al 



FIG. 6; The K vs. Jd phase diagram. Phase boundaries are indicated by dashed hnes and the various trajectories of the 
calculation by solid lines. 




Fig. 7 Path et al. 

FIG. 7: The Oovcn and Oodd string order parameters calculated along path 3. 



